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ABSTRACT 

Aims. We are studying the effects of submm observations on the total dust mass and thus dust-to-gas mass ratio measurements. 
Methods. We gather a wide sample of galaxies that have been observed at subnullmieter (submm) wavelengths to model their Spectral 
Energy Distributions using submm observations (>160 /im) and then without submm observational constraints in order to quantify 
the error on the dust mass when submm data are not available. Our model does not make strong assumptions on the dust temperature 
distribution to precisely avoid submm biaises in the study. Our sample includes 52 galaxies observed at submm wavelengths. Out 
of these, 9 galaxies show an excess in submm which is not accounted for in our fiducial model, most of these galaxies being low- 
metallicity dwarfs. We chose to add an independant very cold dust component (T= lOK) to account for this excess. 
Results. We find that metal-rich galaxies modelled with submm data often show lower dust masses than when modelled without 
submm data. Indeed, these galaxies usually have dust SEDs that peaks at longer wavelengths and require constraints above 160 /xm to 
correctly position the peak and sample the submiUimeter slope of their SEDs and thus correctly cover the dust temperature distribution. 
On the other hand, some metal-poor dwarf galaxies modelled with submm data show higher dust masses than when modelled without 
submm data. Using submm constraints for the dust mass estimates, we find a tightened correlation of the dust-to-gas mass ratio with 
the metallicity of the galaxies. We also often find that when there is a submm excess present, it occurs preferentially in low-metallicity 
galaxies. We analyse the conditions for the presence of this excess and find a relation between the 160/850 //m ratio and the submm 
excess. 
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1. Introduction 



Interstellar Medium (ISM) dust is formed in the ejected ele- 
ments of dying stars. The location of the actual formation of dust 
grains, whether in the ejecta or farther away in the ISM, is still 
debated. Within the ISM, dust undergoes constructive processes 
of accretion or coagulation (Stepnik et al., 2003; Hirashita & 
Omukai, 2009) or destructive/disruptive processes of vaporisa- 
tion, fragmentation, erosion and shattering (Jones et al., 1994, 
1996). Dust grains are, moreover, incorporated in the forma- 
tion of new stars through gravitational collapse of gas and dust 
clouds. Dust grains also participate in the synthesis of molecules 
on their surface (e.g. Hasegawa & Herbst, 1993; Vidali et al., 
2004), leading to the formation of giant molecular clouds within 
the galaxy. Dust finally plays a key role in the overall emission 
of a galaxy since it absorbs stellar and ionized gas radiation to 
re-emit it from Infrared (IR) to submiUimeter (submm) wave- 
lengths, thus shaping its overall Spectral Energy Distibution 
(SED). Dust thus enables us to have access to the most embed- 
ded star fomaing regions of a galaxy which makes it an indirect 
but reliable tracer of its star formation. Dust grains are thus, at 
the same time, relics as well as engines of the evolution of the 
galaxy. 



The ISO and Spitzer telescopes have, in the last few years, 
opened up the view of the interstellar dust in mid-infrared (mid- 
IR) to far-infrared (far-lR) wavelengths. They enabled us to im- 
prove our knowledge of the dust properties of the ISM such 
as the distribution of the different dust components (PolycycUc 
Aromatic Hydrocarbons (R\IIs), silicate, graphite) and how they 
are affected by metallicity (Galliano et al., 2008; Engelbracht 
et al., 2008; Bendo et al., 2008; Munoz-Mateos et al., 2009). 

Many studies have also been performed to investigate the 
relation between the Dust-to-Gas mass ratio (D/G) and metal 
enrichment of galaxies, the D/G of a galaxy being the main out- 
put of dust evolution models and a crucial key to understand 
the cycle of matter within galaxies. Early studies suggested that 
the D/G should be proportional to the metallicity of the galaxy 
(Franco & Cox, 1986). Investigations in the Milky Way (MW) 
found a D/G of ~ 10"^ for our galaxy (Sodroski et al., 1997). 
Further analysis into depletion of metals from the ISM or into 
comparisons of dust extinction versus hydrogen column densi- 
ties have both shown that the D/G of the MW falls within the 
range of 1-5 x 10"^ (see Whittet, 2003; KrUgel, 2003, 2008, 
for more complete details). To investigate the relation between 
the metaUicity and the D/G, Lisenfeld & Ferrara (1998) stud- 
ied a broad sample of dwarf irregular (dirrs) and blue com- 
pact dwarfs (BCDs) from Infrared Astronomical Satellite (IRAS) 
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observations at 60 and 100 yum. A logarithmic correlation be- 
tween D/G and the oxygen abundance was found for the dirrs: 
12 -I- log (O/H) oc (0.52±0.25) X log (MdustP^m), but no rela- 
tion was found for the BCDs which are more actively forming 
stars. Therefore, they observed that metals are less effectively 
incorporated into dust in dwarf galaxies than in more metal- 
rich galaxies. Using the Spitzer Infrared Nearby Galaxies Survey 
sample (SINGS), Engelbracht et al. (2008) also found a steep de- 
crease of the D/G from solar metallicity to a metallicity of 8 and 
note that the D/G seems to flatten out at low metallicity. With 
the same sample, Munoz-Mateos et al. (2009) showed that the 
D/G decreases with galactocentric radius. Bendo et al. (2010) 
also compare radial variations in the D/G with the oxygen abun- 
dance gradients and found a good correspondence, although the 
D/G radial variations seems to be steeper. Recent studies in the 
low-metallicity Magellanic Clouds have shown that their D/G 
are lower than that estimated for more metal-rich environments 
(Bernard et al., 2008; Gordon et al., 2009), consistent with what 
should be expected from the previous relation between the D/G 
and the metalhcity. 

James et al. (2002), using JCMT/SCUBA observations at 
850 fim observations, concluded that the fraction of metals 
incorporated into dust was a universal constant, that is to say 
similar for dwarf and giant galaxies, which was in contradiction 
to the Lisenfeld & Ferrara (1998) conclusions. They also noted 
that their D/G mass ratios estimated with submm observations 
are usually higher than those previously estimated by Lisenfeld 
& Ferrara (1998) for the same galaxies. Galliano et al. (2003, 
2005) also showed that taking the very cold dust traced in 
some low metallicity galaxies by SCUBA or MAMBO, the 
dust-to-metals ratios are closer to the Galactic value. Galametz 
et al. (2009) present the SED models of 4 low-metallicity 
galaxies observed with the APEX/LABOCA instrument at 
870 jum and noticed that the dust masses estimated from their 
models increase when using the submm constraints. On the 
contrary, Draine et al. (2007) found that the dust masses of 
the subset of their SINGS sample observed with SCUBA, 
which are metal-rich galaxies, are rather unaffected by submm 
constraints and that omitting the submm data only moderately 
increases the D/G. Nevertheless, Vlahakis et al. (2005) used the 
SCUBA Local Universe Galaxy Survey (SLUGS) to confirm 
that most of the dust mass of their galaxies was contained in 
cold grains. WiUmer et al. (2009) completed the coverage of 
the far-IR/submm SED of the SLUGS sample with Spitzer data 
and found that these galaxies show flatter vFv(160 yum)/vFy(850 
fim) slopes than the SINGS sample in which galaxies with large 
amounts of cold dust seem to be under-represented. 

Many issues thus still remain open to correctly character- 
ize the D/G with metalhcity, among which is the basic question 
of the accurate quantification of the total dust mass of a galaxy 
and the real impact of submm constraints in the assessment of 
this mass. To compUcate the mass determination issue, there 
have been findings of excess emission in the submm wavelength 
regime (Galliano et al., 2003; Dumke et al., 2004; Galliano et al., 
2005; Bendo et al., 2006a; Marleau et al., 2006; Galametz et al., 
2009; Grossi et al., 2010; O'Halloran et al., 2010; Israel et al., 
2010; Bot et al., 2010, among others). 

Several hypotheses have be studied to explain this excess by 
a variation of the dust emissivity in that submm regime: 

Models of Meny et al. (2007) and observations of Dupac 
et al. (2003) found an effective decrease in the submm emis- 
sivity index as the dust temperature increases. Shetty et al. 
(2009), however, express caution in the inverse temperature - 



/3 interpretation, showing that flux uncertainties, especially in 
the Rayleigh- Jeans regime, can affect the results for the SED 
fits as far as temperature and emissivity are concerned. It is still 
not clear what would be the nature of this new grain population 
or which processes could lead to their predominance in dwarf 
galaxies. Recently, Paradis et al. (2009) interpreted the increase 
of the far-IR emissivity in molecular clouds of the Milky Way 
containing cold dust by fractal aggregates of individual amor- 
phous grains. They note that an increase of a factor of 3-4 in 
the dust emissivity is required to explain the unusually low dust 
temperatures observed by dust aggregation. This explanation of 
a submm excess by grain coagulation was already suggested by 
Bazell & Dwek (1990) (theorethical approach) and Stepnik et al. 
(2001) (using the baloon-borne experiment ProNaos). We note 
that this study is limited to solar metalhcity environments for 
now. Serra Diaz-Cano & Jones (2008) also note that graphite 
can not reproduce the carbonaceous grain sputtering in inter- 
stellar shock waves and cautionned the use of graphite in dust 
models. They proposed that hydrogenated amorphous carbon 
could be the most probable form of carbonaceous grain material. 
Amorphous carbons have a lower emissivity index which could 
be responsible for a flattening of the SED at submm wavelength 
compared to current SED models using graphite (i.e. Meixner 
et al., 2010; Galametz et al., 2010). This scenario alone would 
nevertheless be unable to reproduce the knee observed in the 
submm regime of the SEDs of some low-metallicity galaxies ob- 
served at 450 and 850 jum such as II Zw 40 (see Fig. 3). 

Other explanations try to explain the excess usually observed 
beyond 400 /im by a different population of grains: 

Lisenfeld et al. (2001) and studies of Zhu et al. (2009) on 

extreme extragalactic environments suggested that the submm 
excess could originate from an enhanced very small grain abun- 
dance with an emissivity index /3-1. Ferrara & Dettmar (1994) 
showed that 'spinning dust' located in the ionized gas of many 
galaxies were producing radio emission. Studies by Anderson 
& Watson (1993), Draine & Lazarian (1998), Hoang et al. 
(2010), Silsbee et al. (2010) or the recent Planck results (Planck 
Collaboration, 2011a) aimed to characterise this component 
which explains most of the 10-100 GHz emission, also called 
the anomalous foreground. Ysard & Verstraete (2010) recently 
suggested that their peak frequency could depend on different 
parameters (radiation field intensity, size distribution of the dust 
species, dipole moment distribution, physical parameters of the 
gas phases etc.) and thus that spiiming dust emission could be re- 
sponsible for the excess at submm wavelengths observed in some 
galaxies, as recently suggested by Bot et al. (2010) for the Large 
Magellanic Cloud. Spinning dust emission being generated by 
very rapidly spinning and assymetric small dust grains, PAHs 
are often proposed as possible carriers. However, this assump- 
tion seems to be in contradiction with the lack of PAH features 
observed in low-metallicity galaxies in which submm excess is 
usually detected. 

This excess could finally be interpreted as a very cold 
(<10K) dust component (e.g. Galhano et al., 2003, 2005; 
Galametz et al., 2009). This latter hypothesis could explain the 
break observed in the submm regime of some SEDs, especially 
in low-metallicity galaxies. However, its use in SED models 
usually implies a significant increase of the total dust mass, with 
more than 50% of the total dust mass residing in this very cold 
component, leading to unrealistic D/G in some galaxies (e.g. 
Dumke et al., 2004; Bendo et al., 2006b), if we are taking into 
accoimt the total gas reservoir detected. 
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The ISO and Spitzer telescopes offered a vision of the IR 
emission limited to 160 /zm and thus do not place constraints 
on the bulk of the coldest (<15K) dust of the galaxies, that 
could represent a large amount of the dust mass of galaxies. For 
many galaxies, the Herschel telescope, more specifically the 
SPIRE instrument observing from 250 to 500 yum, will better 
constrain the subnun slope and thus help to disentangle between 
those scenarios. The issues discussed here will continue to be 
of concern throughout the Herschel, Planck and ALMA era, 
especially for low-metallicity objects not detected. Thus, this 
study lays the ground work for an eventual broader study from 
which these future data can build on. 

We perform a systematic study of the dust SEDs of galax- 
ies spanning of a wide range of metalhcities with the following 
objectives: 

1) compare the estimates of the total dust mass of the galax- 
ies when the SED modelling has submm constraints, to the dust 
masses determined when the subnmi observations are not used 
(as if they are not available). 

2) using only the sample modeled with submm data, to have 
a reliable dust mass, determine the D/G as a function of metal- 
licity and compare with chemical evolution models. 

3) identify the galaxies of the sample which show a submm 
excess. 

This paper is organised as follows: in § 2, we describe the 
preliminary study we perform on the D/G as a function of metal- 
licity by gathering dust and gas masses from the hterature for a 
large sample of galaxies. In § 3, we describe the coiimion SED 
modelling technique we apply to remodel a subsample of these 
galaxies in a unique way with and without submm constraints. 
Finally, in § 4, we analyse the differences obtained in the dust 
masses found with or without submm constraints and analyse 
the conditions under which a subnnm excess is found in galaxies 
and relate this to the metallicity. We finally summarize the main 
conclusions of the study in § 5. 

2. Genesis of our study 

For this §, we perform a preliminary study based on data directly 
taken from the literature, in particular, the dust masses. 

2.1. The sample 

We have gathered published dust and gas masses of a broad sam- 
ple of galaxies, some of them observed with submm telescopes, 
to study the relation of D/G with metallicity. The dust and gas 
masses used in this preliminary study are published in: 

Lisenfeld & Ferrara (1998): This sample groups together dwarf 
irregular galaxies (dirrs) and BCDs. Their SEDs were per- 
formed using only IRAS data and a single temperature, mod- 
ified black body model. Their dust mass estimates are thus 
considered as lower hmits to the total dust masses of the 
galaxies. 

James et al. (2002): In this sample, the submm observations 
were performed with JCMT/SCUBA at 850 jim. They fitted 
their data with the Dunne & Eales (2001) two-component 
modified blackbodies to derive the total dust mass of their 
galaxies. 

Hirashita et al. (2008): This sample of eight low-metallicity 
galaxies was observed with AKARI at 65, 90, 140 and 160 
jum. They derived a total mass of dust for their galaxies 



from the far-IR emission using a single temperature modi- 
fied blackbody model, using the mass absorption coefficient 
of dust grains taken from Hildebrand (1983). 

Draine et al. (2007): Their sample (SINGS, see Kennicutt et al., 
2003) includes spirals, ellipticals, starburst galaxies and both 
metal-rich and metal-poor galaxies. Their dust SED mod- 
elUng use the Li & Draine (2001) dust properties and is per- 
formed with Spitzer IRAC and MIPS data, with some of their 
galaxies observed with SCUBA at 450 or/and 850 jim. 

Engelbracht et al. (2008): Their galaxies, observed with 
Spitzer, cover a wide range of metallicity values (from 
12-i-log(0/H)=7.1 for SBS 0335-052 to 8.85 for IC342). 
Their dust masses are estimated using the standard formula 
and the absorption coefficients from Li & Draine (2001). 

Galliano et al. (2008): A subset of their sample had been ob- 
served with SCUBA 450 /zm and/or 850 fan. All of their 
SED models included either the SpitzerflRS spectra or the 
ISOCAM spectra for optimum contraints in the mid-IR. 
Their model includes a detailed decomposition of dust emis- 
sion into its gas-phase components. 

Galametz et al. (2009): their sample contains 3 low-metallhcity 
sources observed with the Spitzer telescope and LABOCA, 
a 870 yum multi-channel bolometer array available on the 
APEX telescope in Chile. 

In these papers, the gas masses are estimated using Hi 21cm 
line measurements and the H2 mass determined from CO obser- 
vations, when available. 

2.2. The relation between the dust-to-gas ratio and 

metallicity 

Figure 1 shows the D/G of these galaxies as a function of metal- 
licity. We overlay several chemical evolution models on this 
figure: the one-zone single-phase chemical evolution model of 
Dwek (1998) presented in Galliano et al. (2008) and the dust 
formation models of Edmunds (2001). Edmunds (2001) note 
that they adopt a less detailed quasi-analytic model compared to 
Dwek (1998). Their models include mantle growth in the ISM 
and take into account the grain core production in both super- 
novae and giant stars. These two modes of dust formation are 
plotted in Fig. 1 . The solid line corresponds to dust produced in 
low and intermediate mass stars and the dotted fine corresponds 
to dust production strongly affected by the role of supemovae. 

The filled symbols localise the galaxies for which dust 
masses are estimated using submm constraints. These galaxies 
tend to follow the Edmunds (2001) and Dwek (1998) models 
and clearly differ from the linear regression performed with the 
whole sample (dashed line), i.e. which also include galaxies 
that were not observed at submm wavelengths. We particularly 
noticed that low-metallicity galaxies (12-i-log(0/H)<8.5) whose 
dust masses were estimated with submm constraints (filled sym- 
bols) have higher D/G than galaxies modelled without submm 
constraints (empty symbols, crosses or asterisks). This could 
imply that for dwarf galaxies, the use of submm constraints 
leads to higher dust mass predictions. We further investigate this 
behaviour in this study. 



2.3. The limitations of the samples 

Of course, these various samples do not necessarily have any 
selection criteria in common. Some of the surveys were de- 
liberately targeting dwarf galaxies yet others are covering a 
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wider variety of galaxies. For instance, the SINGS sample is 
designed with a broad range of galaxies evenly filUng a three- 
parameter space defined by morphology, infrared luminosity and 
LfiR/Lopticai- The dust models and SED fitting procedure used to 
derive D/G are different from one sample to another, which could 
lead to biases in the analysis and interpretation. Moreover, all of 
the galaxies used in this broad sample have not been observed 
at submm wavelengths, data that are necessary to put constraints 
on the coldest phases of the dust which could represent a large 
amount of the total dust mass. Such a global study should quan- 
tify the total gas mass in a uniform manner. We discuss this issue 
in Section 4.2.1. 

In conclusion, analyzing the dust and gas masses of these 
galaxies directly given by the literature as one broad collection 
should be limited. We decided 1) to reprocess this analysis us- 
ing a common realistic SED model to avoid biases in comparing 
galaxies analysed with different SED models; 2j to restrict our 
study to galaxies observed at submm wavelengths to determine, 
for a set of fixed dust properties (emissivity, composition), the 
influence of submm constraints on the estimates of the total dust 
mass of the galaxies. Thus we use a subsample of the galaxies 
from Fig. 1 in the following study, namely galaxies which have 
submm observations. 

3. The SED models 

3.1. The wavelength coverage 

From the previous full sample described in § 2, we only keep 
galaxies for which submm data are available, namely observa- 
tional constraints beyond 200 //m. The restricted sample is pre- 
sented in Table 1 and gathers subsamples of the James et al. 
(2002), Draine et al. (2007), Galliano et al. (2008) and Galametz 
et al. (2009) samples. 

MetaUicities, distances and submm 450 and/or 850 (or 870) 
jum observations for this subsample are summarized in Table 1. 

3.1 .1 . Submm flux densities 

Most of the submm data we used are flux densities at 450 
and/or 850 yum from SCUBA or 870 yum from LABOCA. A 
few galaxies of our sample were not observed with SCUBA 
or LABOCA but with other submm instruments. M83 was ob- 
served at 540 fim with the 4-m telescope at Cerro Tololo Inter 
American Observatory (Hildebrand et al., 1977). NGC 1068 was 
observed at 350 fim and 390 /im, respectively with the Submm 
High Angular Resolution Camera (SHARC) (Benford, 1999) 
and with the University of Hawaii 2.2-m telescope at the Mauna 
Kea Observatory (Hildebrand et al., 1977). NGC 891 was ob- 
served at 350 and 1300 fim, with the 3-m Infrared Telescope 
FaciUty (IRTF) and the 88" telescope of the University of Hawaii 
on Mauna Kea (Chini et al., 1986). Finally, IC 342 was observed 
at 1mm with the Palomar 5-m telescope (Elias et al., 1978). Only 
an upper limit was derived for this galaxy. 

We would also like to wam the reader about the quality 
of the subtmn fluxes derived for some of the SINGS galaxies 
used in this paper. Indeed, Draine et al. (2007) warned that the 
galaxy NGC 6946 was scan mapped with SCUBA and that the 
data processing could have removed some extended emission 
for this galaxy, leading to an underestimate of submm fluxes. 
Furthermore, in the case of NGC 1097, SCUBA observations 
only cover the galaxy center and the total submm fluxes were es- 
timated by inferring how much flux resides in the outer regions 
not observed, assuming a constant MlPS-to-submm ratio (Dale 



et al., 2007). This assumption could be incorrect if radial color 
variations are present. SCUBA maps of NGC 3521, NGC 3627, 
NGC 4536 and NGC 7331 also have small field of views. The 
submm fluxes presented in Dale et al. (2007) should therefore be 
used with caution in a global analysis of the galaxies. 



3.1 .2. Optical and infrared flux densities 

2MASS data constrain the stellar contribution of the SEDs. 
2MASS flux densities are taken from the 2MASS Large Galaxy 
Atlas (Jarrett et al., 2003), Galametz et al. (2009) and Dale et al. 
(2007). Spitzer/IRAC and MIPS flux densities are used to sam- 
ple the infrared part of the SED, when available. The Spitzer flux 
densities used in the following study are also given in Dale et al. 
(2007), Engelbracht et al. (2008) and Galametz et al. (2009). 
For the SINGS galaxies used in this paper and for which far-IR 
data are critical when submm data are not available, the Spitzer 
data enable us to sample both the peak position and the begin- 
ning of the submm slope. For the Galliano et al. (2008) sample, 
SpitzerflRS spectral information is, in most cases, available and 
used in our modeUing. These spectra are presented in Galliano 
et al. (2008). 

We use file best available IRAS fluxes. The IRAS flux den- 
sities of the Draine et al. (2007) galaxies were obtained from 
SCANPI and the IRAS High Resolution Image Restoration 
Atlas (HIRES; Surace et al., 2004). Those of the James et al. 
(2002) sample are coming from the IRAS Revised Bright 
Galaxy Sample (Sanders et al., 2003), HIRES for the galax- 
ies NGC 3994 and NGC 3995 and from flie IRAS Faint Source 
Catalog (Moshir & et al., 1990) for NGC 4670. The IRAS flux 
densities of Galhano et al. (2008) and Galametz et al. (2009) 
galaxies are finally obtained from the IRAS Faint Source cata- 
logue and the IRAS catalog of Point Sources (IPAC 1986). The 
choice of the IRAS data is of little influence on the dust masses 
we derive from the SED models. We finally note that there is no 
systematic difference between the wavelength coverage of dwarf 
galaxies compared to higher metaUicity galaxies. 



3.2. The building blocks of our fiducial SED model 

The SED of a galaxy synthesizes the emission at each wave- 
length of its difl'erent components: H n regions, molecular 
clouds, diflrise ISM etc. The SED fits obtained from any model 
are dependant on the assumptions we make a priori about pa- 
rameters such as radiation field or dust properties. Despite this 
dependance, valuable comparisons can be made from the study 
of different galaxies using the same SED modelling method. 

We first consider that the sources of IR emission are the dust 

grains (PAHs, silicate, graphite) and the stellar continuum. Our 
SED modelling approach uses realistic intrinsic dust properties 
and is commonly used in the litterature (Dale & Helou, 2002; 
Paradis et al., 2009; Whaley et al., 2009). The properties 
of silicate and grapite grains, namely their size distribution 
and chemical composition, are the standard and up-to-date 
properties of the Bare-Gr-S dust model Zubko et al. (2004). We 
emphasis that this model is consistent with the recent results of 
Draine et al. (2007) and the DustEM model of Compiegne et al. 
(201 1). We assume that the dust particles are thus a mixture of 
grains with solar depletion constraints. The optical properties 
of the sihcate and graphite grains are taken from Weingartner 
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Fig. 1. Dust-to-gas mass ratio as a function of oxygen abundance. Black crosses show the Lisenfeld & Ferrara (1998) sample of Blue Compact 
Galaxies and dwarf Irregulars. Purple asterisks show the Engelbracht et al. (2008) sample using Spitzer data. Blue squares indicate the James et al. 
(2002) galaxies. Green circles show the galaxies of the SINGS survey (Draine et al., 2007) for which metallicity is already published. Orange 
upward triangles are the Galliano et al. (2008) sample, blue downward triangles are the Hirashita et al. (2008) sample. Stars finally indicate the 
location of the three low-metallicity galaxies Haro 11, NGC 1705 and Mrk 1089 of Galametz et al. (2009). Filled triangles, squares, circles or stars 
indicate when submm data were used in the SED modelling to estimate the total dust mass of the SED. The solid line and the dotted line are the 
predictions of the dust formation models from Edmunds (2001): the solid line shows a model in which dust is produced in low and intermediate 
mass stars while the dotted line shows the model in which supernovae play an important role in the production of dust. The gray stripe is the 
expectation of the one-zone, single-phase chemical evolution model of Dwek (1998) presented in Galliano et al. (2008). The dashed line shows 
the linear regression of the whole sample. The D/G errors can be significant for many galaxies but should not exceed ± 50%. (the ± 50% error bars 
are plotted on the bottom right part of the plot) 



& Draine (2001) and Laor & Draine (1993) respectively \ 
In order to better constrain the SED in the mid-IR range, we 
introduce two free parameters to model the PAH emission: the 
PAH charge fraction ifpAH+) and the PAH-to-dust mass ratio 
(fpAH), normalised to the Galactic value of 0.046. The PAH 
and PAH-H properties (e.g. the absorption efficiencies) are taken 
from Draine & Li (2007). When observations are lacking 8 fj.m 
measures and to avoid a degeneracy in the modelling, this fp^H 
parameter is fixed to 0.5 (which means half the Galactic value). 
This choice does not significantly influence the total dust mass 
derived from flie models. The total mass of dust Mjust is a free 
parameter of our modelling. 

We adopt the spectral shape of the Galactic diffuse ISM 
Interstellar Radiation Field (ISRF) determined by Mathis et al. 
(1983). Previous studies have shown that the form of the 
global ISRF is different from normal or dusty galaxies to low- 
metallicity galaxies (e.g. Madden et al., 2006), an effect at- 
tributed to the lower dust attenuation and thus the larger mean 
free path length of the ionizing photons in low-metallicity ISM 
and the predominance of massive stars in those environments. 
Nevertheless, the shape of the ISRF does not influence the far- 
IR to submm regime of the SED and thus the total dust mass 
derived (see Galametz et al., 2009). The radiation field indeed 



' We also refer to Draine & Li (2007) and 
http://www.astro.princeton.edu/~draine/dust/dust.html where dust 
properties can be found. 



mainly affects the emission of the non-equilibrium grains that 
do not significantly contribute to the dust mass. 

The SED is multiplied by a dimensionless scaling factor U, 
with U-l corresponding to a radiation field of the local solar 
neighborhood value of 2.17 x lO"'* W m"^ (Mathis et al., 1983). 
To describe the various environments within which the various 
components of the ISM reside, we combine the individual 
SEDs using the prescription defined in Dale et al. (2001). 
This parametrisation of the analytical form of the temperature 
distribution has been widely adopted to describe the ISRF 
in galaxies, for example by Dale & Helou (2002), Galliano 
et al. (2008), Paradis et al. (2009), among others. It assumes a 
power-law distribution of the dust mass over a wide range of 
interstellar radiation heating intensities and conservatively links 
the hot and cold emission in galaxies: 



dMdustiU) cc U-"dU (1) 

where M^j,,, is the dust mass heated by the radiation field, U its 
intensity and a the coefficient of the power law describing the 
contributions of the individual SEDs (also a free parameter of 
our model). 

Radiation field intensities are ranging between the two free 
parameters U,„,„ and {]„,„,[. The behavior of the emission of 
the different dust components as a function of U are plotted 
in Fig. 2. We note that Draine et al. (2007) used a similar 
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Fig. 2. Emission of silicate, graphite, neutral and ionized PAHs for different values of U. 
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parametrisation but add a "diffuse ISM" component with a 
unique intensity factor \J-\J mm- They choose to restrain this 
parameter to values > 0.7 as an ad hoc measure to avoid 
invoking large dust masses when the submm slope of the SED 
is not constrained. This assumption is further discussed in § 4.1. 

Finally, we subtract the stellar contribution to the SED in the 
mid-IR using 2MASS and the SpitzerflRAC 3.6 fim bands. We 
model the contribution of old stars using the stellar evolution 
program PEGASE (Fioc & Rocca-Volmerange, 1997), assuming 
that the stellar population has undergone an instantaneous burst 
5 Gyr ago. We adopt the Salpeter Initial Mass Function : N(M) 
dM oc M"' ■'^ dM (Salpeter, 1955), and an initial metallicity 
of Z-Zq. The mass of old stars Moidsiars is a free parameter 
of our modelling. The choice of the IMF does not particularly 
influence our study since we are only fitting the IR and do not 
want to necessarily quantify the mass of old stars. We finally 
note that no radio component was introduced in the modelling 
compared to the similar model of Galliano et al. (2008). 



The overall modelling is an iterative process. We compute 
the temperature distribution of the dust grains heated by the ab- 
sorption of the stellar radiation. The IR emission is then precal- 
culated for each dust grain (different sizes, different properties) 
and each intensity U. The sum of these discrete contributions 
leads to a global SED model of the galaxy: 



with 

j^PAH+diisI 



■'(U,A)U-''M,,.„dU + LT''{m) 



+ Lf'"'''''"{U, A) + Lf"""{U, A) 



(3) 



In summary, the free parameters of our modelling are: 

- Mj„s, the total mass of dust 

- fpAH the PAHs-to-dust mass ratio (normalised to the Galactic 
value) 

- fpAH+ the ionized PAHs-to-total PAHs mass ratio 

- a the index describing the fraction of dust exposed to a given 
intensity 

- U,„,„ the minimum heating intensity 

- UmoA the maximum heating intensity 

- Goldstar the mass of old stars 



Both IRAS and Spitzer bands are color-corrected. The IRAS 
fluxes are corrected using a v x Fy convention. The Spitzer 
bands are corrected using the convention described on the 
Sp/teer website ^. 

Finally, all data are weighted with respect to the density of 
their neighbours. The model obtained is integrated within each 
filter and compared to observed values using a minimisation 



^ IRAC: http://ssc.spitzer.caltech.edU/irac/iracinstrumenthandbook/2 1/ 
MIPS: http://ssc.spitzer.caltech.edu/mips/mipsinstrumenthandbook/51/ 
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based on the Levenberg-Marquardt algorithm (least squares fit- 
ting routines): 



X 



n- 



2/1; 



Lf^Ad - LMi) 



ALfHAi)l2 



(4) 



where A L,f'{Ai) is the error on the luminosity at a given 
wavelength. 



3.3. Modelling the submm excess 

Most of the galaxies of our sample can be fit with our SED 
model up to submm wavelengths. Nevertheless, for 9 galaxies 
(namely NGC 1569, NGC 1140, He 2-10, II Zw 40, NGC 4303, 
NGC 2903, NGC 7674, NGC 1705 and Haro 11), our model 
can not fit together the far-IR measures and the elevated submm 
fluxes and an excess at submm wavelengths is detected. Those 
galaxies are mostly low-metallicity dwarf galaxies. Such excess 
is often reported while modelling the SEDs of low-metallicity 
galaxies (Galliano et al., 2003; Dumke et al., 2004; Galliano 
et al., 2005; Bendo et al., 2006a; Marleau et al., 2006; Galametz 
et al., 2009; Grossi et al., 2010; O'Halloran et al., 2010; Israel 
et al., 2010; Bot et al., 2010). 

We note that Li & Draine (2001) modified the standard siU- 
cate emissivity at A> 250 fim in order to match the average high 
Galactic latitude dust emission and the Galactic "submm excess" 
observed with DIRBE/FIRAS. In this paper, we use the Zubko 
et al. (2004) model that fits the same constraints. The submm 
excess detected in our 9 galaxies may have the same origin as 
the Galactic submm excess, but the amplitude of the excess ob- 
served in the Milky Way is much smaller than observed in dwarf 
galaxies. 

This paper will assume that the excess is due to cold dust 
and explore the consequences of this hypothesis on the global 
properties (total dust mass and D/G) of those 9 galaxies. For 
these 9 galaxies, we thus add an independent thermal cold dust 
component to our fiducial model. The luminosity of the cold dust 
component is calculated using: 



/A X-^ 

L^^(i) = McDKy,cD [-jj (47tBAA,Tcd)) 



= McD 



(4/3 ..3,) (y) (4-^v(^^^c.)) 



(5) 



where By is the Planck function. Tea the cold dust temperature,/? 
the emissivity coefficient. We adopt Aq-\QQ fim for the reference 
wavelength. We assume graphitic values for the mass density of 
the grains p=2.5 x lO-' kg m"-' (close to the mass density of 
carbonaceous grains of 2.2 x lO-' kg m"^ of Draine & Li, 2007) 
and the absorption coeSicient Qo/a =150 cm ' at Aq, with a, the 
radius of the grain. We are studying an effect which is difficult 
to constrain completely with one or two submm fluxes. 

The lack of constraints forces us to fix the temperature 
and the dust emissivity index of this blackbody to conservative 
values. 

Temperature - Galliano et al. (2003, 2005) and Marleau 
et al. (2006) suggest that their submm excess observed in 
low-metallicity galaxies could be explained by the presence of 
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Fig. 3. SED models of the galaxies for wliicli an excess at submm wave- 
lengthis is detected. The SEDs modelled with submm data are indicated 
by the black lines while the dashed black lines indicates the model ob- 
tained when submm data are omitted in the modelling. Observational 
constraints are superimposed (filled circles). When the error bars are 
not shown, the errors are smaller than symbols. The thick orange line 
shows the IRS spectrum used in the SED modelling when available. 
The open circles indicate the expected modelled fluxes integrated over 
the instrumental bands. The green, red and blue lines respectively dis- 
tinguish the stellar, the warm dust and the cold dust contributions to the 
SED determined with submm data. 



cold dust distributed in small and dense clumps and modelled 
the excess with an independant cold dust component with 
temperatures between 5 and 10 K. We tried modeled the 9 
galaxies with a cold dust component with T>10K (ex:15K). 
In most cases, this leads to fits with higher chi-square values. 
When the 450 fim measure is available, the model has difficulty 
to fit the peak of the SED and the 450 fim data together and 
usually leads to a SED model that overestimates the 450 //m 
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Fig. 3. continued 



data. The temperature choice influences the dust mass derived 
from our modefling and we decide to fix the cold dust compo- 
nent temperature to lOK. Using lOK (and not <10K) dust also 
prevents us from requiring excessive cold dust. 

Emissivity - A flattening of the submm slope is often ob- 
served when 450 + 850 or 870 //m measurements are available 
(see Fig. 3). As suggested by several hypothesis to explain the 
submm excess and reviewed in the introduction of this paper, an 
increase of the emissivity at submm wavelengths is expected to 
explain the flattening of the submm slope. Moreover, Galliano 
et al. (2003, 2005) and Galametz et al. (2009) modeled the 
submm excess they detected with a very cold dust component 
and found that an emissivity index of 2, as commonly used to 
describe the submm slope of SEDs (i.e. Dunne & Bales, 2001; 
Draine & Li, 2007) leads to unrealistically high dust masses. 
(For instance, Galametz et al., 2009, found that dust masses 



were multiplied by a factor of 6-7 when using [3-2 in lieu of 1). 
We conservatively adopt the value fi-l. We fix our absorption 
to the model of Li & Draine (2001), which may deviate when 
a dust emissivity index of 1 (and not the standard emissivity of 
2) is used. We assume that this value would still be applicable 
here. 
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The SEDs of the 9 galaxies showing a submm excess and 
modelled with an extra cold dust component are presented in 
Fig. 3. We already want to note that those galaxies only consti- 
tute a small subset of our sample and global relations studied 
in the following should not be affected by those objects. Those 
galaxies are marked by red symbols in the following plots. 



4. Analysis 

4.1. Dust mass estimates with submm observations 

We apply our SED modelling technique to the observational con- 
straints we gathered. We derive dust mass estimates with all the 
data including the submm data, and another time omitting the 
submm data, which normally means that the observations stop 
at 160 fim. To quantify the errors on our dust masses for each 
SED solution, we produce a grid of 100 randomly modified ob- 
servational constraints, allowing the observational fluxes to vary 
within their error bars, following a Gaussian distribution around 
their reference value. Dust masses and errors are summarized in 
the last column of Table 1. 

We also test the presence of the mid-IR spectroscopic ob- 
servations (e.g. ISO or Spitzer IRS) on the dust mass deduced 
without submm, to study the influence of submm data on the 
submm slope. We run our model with and without the mid-IR 
spectrum and without data beyond 200 /im. We find no strong 
differences nor systematic increase or decrease of the total dust 
mass deduced when a mid-IR spectrum is added to constrain the 
SED, as long as there are at least a few observations at mid-IR 
wavelengths (i.e. IRAC and MIPS 24 iim or IRAS 25 fim). 

In Fig. 4a, we show the values of the dust mass estimated 
with submm data vs the dust mass estimated without submm 
data. Figure 4b presents the ratio of these two quantities as a 
function of metallicity. Both figures demonstrate that, for dusty 
massive galaxies, the dust mass tends to be higher when submm 
measures are omitted. We also note that errors on the dust mass 
of galaxies (as given in Table 1) can be significant when submm 
data are omitted, especially for dustier galaxies. 

In the case when submm data is lacking, there is indeed a 
risk that the SED procedure invokes a large mass of cool dust. 
Fig. 5 illustrates this effect on the galaxy NGC 337 modelled 
with (solid line) and without (dashed line) the 850 yum SCUBA 
flux. The submm slope of the SED is reasonably sampled when 
the submm flux is available and the total dust mass decreases 
by a factor of 8 for this galaxy. This systematic overestimate 
of the dust masses mostly affects metal-rich galaxies. Indeed, 
their SEDs usually peak at longer wavelengths compared to 
the low-metallicity galaxies. The 160 flux is thus no longer 
sufficient to fully sample the dust SED peak and constrain the 
submm slope. 

The SCUBA/LABOCA fluxes constrain the submm slope of 
our SEDs and more particularly the minimum intensity Um,„. 
This parameter is directly linked to the cold dust mass (see 
Fig. 2 and the influence of <U> in the emission of silicate 
and graphite) and thus significantly influences the total dust 
mass derived. An overestimate of U,,,,-,, will systematically 
lead to an underestimate of the dust mass. Figure 6 shows the 
evolution of U,„,>, parameter if submm observations are used 
in the modelling. The galaxies showing a submm excess do 
not appear on this plot. We observe a tighter dispersion of this 
parameter when submm data are used in the modelling. The 
small values of this parameter when submm observations are 
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Fig. 4. a) Total dust mass estimated witli submm observations versus 
total mass of dust estimated without submm observations. The line 
shows the 1 : 1 relation. Squares indicate the James et al. (2002) galaxies. 
Circles indicate the galaxies of the SINGS survey (Draine et al., 2007) 
for which metallicity is already published. Triangles are the Galliano 
et al. (2008) sample and stars indicate the location of the three galaxies 
Haro 11, NGC 1705 and Mrk 1089 of Galametz et al. (2009). Red sym- 
bols show the galaxies that present a submm excess and for which a cold 
dust (lOK) component with an emissivity /3=1 is added, b) Evolution 
of the ratio of the total mass of dust estimated with and without submm 
data as a function of the oxygen abundance of the galaxy. 



omitted lead to the high dust mass we derive in some cases. 
A similar study on the Um,„ parameter has been performed 
by Draine et al. (2007). They found that the dust masses of 
their SINGS-SCUBA sample is rather unaffected by the use 
of submm constraints. For the same galaxies, we find, on the 
contrary, that the dust masses we derive without using submm 
data can be significantly higher than those obtained using 
submm data. Our model assumes a power-law distribution of 
the dust mass over interstellar radiation heating intensities. 
The SED model used in Draine et al. (2007) includes both 
this parametrisation and a "diffuse" dust component heated 
by a single radiation field with an intensity of U=U,„,„ which 
describes most of the far-IR / submm emission of their galaxies. 
They also chose to restrict their radiation fields to intensities 
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Fig. 5. Total SED of the galaxy NGC 337 using the 850 /^m SCUBA 
flux (solid line) and not using the submm information (dashed line). 
Observational constraints are superimposed (filled circles). Open cir- 
cles are the expected modelled fluxes integrated over the instrumental 
bands. Green and red lines respectively distinguish the stellar and the 
dust contributions of the SED. 



100.0 



S 1.0 








1 ' 


n 





r 










. / 




















• • 










• 


















t 








_ _ _ 










t 


• 














• < 








1 


< >' 1 




1 





0.001 0.010 0.100 1.000 10.000 lOO.OOC 
Umin without submm constraints 

Fig. 6. Minimum heating intensity derived from the SED model using 
submm data compared to minimum heating intensity determined with- 
out submm data. The orange line indicates the 1 : 1 relation. The dashed 
lines define the limits Draine et al. (2007) used as boundaries for their 
minimum heating intensity U„,i„. We remind the reader that our model 
solutions are obtained with U^m as a free parameter. For convention 
on symbols, see Fig. 4. Limiting U„„„ to 0.7 avoids an overestimate of 
the dust mass but removes the possibility of low temperature radiation 
fields. 

U„„„ > 0.7 as an ad hoc measure to avoid their model invoking 
a large mass of cold dust heated by weak starlight when submm 
observations are not available. The boundaries of the minimum 
heating intensity range chosen by Draine et al. (2007) are shown 
in the shaded region in Figure 6. The restriction of the U„„„ 
interval seems to be consistent with their objects, and for most 
of the metal-rich galaxies of our sample. 

However, as noted in the Draine et al. (2007) study, the 
SINGS-SCUBA galaxy sample could be biased in favor of in- 
creased dust mass and associated star formation. The restriction 
on the Um,„ parameter would not be valid in systems containing 
a large amount of cold dust heated by weaker stellar light. This 



assumption could lead, in a systematic way, to an underestimate 
of the dust mass of those environements when submm data are 
not available and should be used with caution in some galaxies. 

Indeed, we find that for some low-metallicity galaxies of our 
sample, the dust masses estimated without submm observations 
are, on the contrary, lower than those determined with submm 
data (Fig. 4b). This increase of the total dust mass while us- 
ing submm data was already discussed in James et al. (2002). 
They observed some of the galaxies presented in Lisenfeld et al. 
(2001) with SCUBA and noticed an increase of the total dust 
mass of their galaxies when submm data are used. Indeed, the 
Lisenfeld et al. (2001) study was performed based on IRAS 
60/100 /vm fluxes only and could have missed cold dust at longer 
wavelengths and thus underestimated the total dust mass. For 
most of the low-metallicity galaxies of our sample, the restric- 
tion of the U„„„ parameter when submm data are not available is 
an unsatisfactory solution, since such a bias leads to an underes- 
timate of the total dust mass. In those environments, submm data 
are crucial to trace the coldest phases of dust and also detect, if 
present, a possible submm excess. 

In conclusion, when comparing the results unconstrained to 
those with the proper submm data, we find that the uncertainty 
on the dust mass is non-trivial and systematic, depending on 
metallicity of the object. Without submm constraints, our fit of 
the SEDs will give higher mass for metal-rich galaxies and lower 
masses for low-metallicity galaxies. 

In Fig. 4, we finally note that our modelling of the submm 
excess (which includes a modified blackbody with T=10K - red 
symbols) does not systematically lead to an increase of the dust 
mass estimate of a galaxy compared to that estimated omitting 
submm data . The slope of the trend (i.e. the correlation coeffi- 
cients) observed in Figure 4b is, furthermore, not affected by the 
addition of galaxies in which a submm excess is detected. 

4.2. The Dust-to-Gas mass ratio relation witti metallicity 

4.2.1. The gas masses and metallicities 

Values and references of the gas masses of our galaxies are 
given in Table 1. Both atomic and molecular gas masses are 
given when available in the literature and gas masses were eval- 
uated using the same distances as those used to derive the dust 
masses and quoted in Table 1. We determine the molecular gas 
masses using the canonical galactic Xco factor of 2.3 x 10^" 
cm-2 / K km s"' of Strong et al. (1988) (Table 1, Value g«/). The 
conversion factor from CO to H2 (Xco) in galaxies depends on 
the metallicity (Wilson, 1995; Israel, 1997, 2000; Barone et al., 
2000; Boselli et al., 2002; Strong et al., 2004; Leroy et al., 2005). 
Boselli et al. (2002) studied Xco in normal late-type galaxies and 
derived Xco factors ranging from 10^° cm"^/K km s"' in giant 
spirals to 10^' cm"^/K km s"' in dwarf irregulars. They conclude 
that using a constant Xco factor to derive the molecular gas mass 
from CO observations could lead to an underestimation of a fac- 
tor ~10. We use the relation from Boselli et al. (2002) to estimate 
Xco for each galaxy depending on its metallicity and rescale the 
molecular gas masses (Table 1, Value ^02): 

logXco = -1.01(12 -I- Iog(0/H)) + 29.28 (6) 

The metallicities of the galaxies are given in Table 1 and 
were estimated in the different papers using standard R23 meth- 
ods ^, except in some cases discussed in § 4.2.2. We would like 

^ The oxygen abundance, commonly used to estimate the metallic- 
ity, can be derived from the bright line observations and depends on 
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to warn the reader about the current debate on the correct way 
to estimate the metallicity of a galaxy. For instance, Pilyugin & 
Thuan (2005) discuss the problem of empirical calibrations per- 
formed with strong oxygen line intensities-oxygen abundances 
and propose revised methods to estimate the metallicity (see 
Moustakas et al., 2010, for applications on the Spitzer/SINGS 
sample). These different methods will not affect the general trend 
of our sample but could, nevertheless, introduce a systematic off- 
set in our metallicity values. 

4.2.2. A revised D/G relation 

Figure 7 shows the D/G as a function of oxygen abundance. 
More precisely, the D/G is estimated using the dust masses cal- 
culated without subirnn fluxes in Fig. 7a and then recalculated 
using the submm data in Fig. 7b. The error bars indicate the 
errors on the dust masses. We separate the two plots for clarity. 
Figure 7c gathers information from Figure 7a and Figure 7b. 
The dot-dash vertical lines indicate how galaxies "move" on 
the D/G axis when subirnn measures are used in the dust mass 
estimation and clearly show that for dustier galaxies, submm 
data significantly affect the total dust mass estimates by shifting 
the dust masses to lower values. The predictions of the Edmunds 
(2001) and Dwek (1998) models presented in § 2 are indicated 
on the plots. We also overlay the linear regression (dashed 
line) which was determined from our original broad sample of 
galaxies (§ 2 and Figure 1). 

Our sample restricted to galaxies observed at subimn wave- 
length does not follow the regression of the original big sam- 
ple presented in Figure 1. The D/G vs metallicity relation, in 
fact, now tends to flatten toward lower metallicity when submm 
data are used in the dust mass estimates, with D/G values of 
low-metallicity galaxies systematically higher using submm data 
than those determined without submm data. Furthermore, the use 
of submm data significantly tightens the D/G vs metallicity re- 
lation. In order to quantify the scatter in both plots, we calculate 
the standard deviation of log(D/G) in the three metallicity bins 
Z=[8.2 , 8.5], Z=[8.5 , 8.8] and Z=[8.8 , 9.1]. The standard devi- 
ation is equal to 0.55, 0.44 and 0.72 in those different metallicity 
bins when submm data are omitted in the SED modelling to de- 
termine the total dust mass. Those values decrease to 0.34, 0.32 
and 0.5 if submm data are used, consistent with the visual im- 
pression of a tightening in the D/G vs metallicity relation. 

We note that galaxies modeled with an extra cold dust 
component (red symbols) do not affect the global relation. 
Unfortunately, the current observations are still not sufficient 
to discriminate between the two models of Edmunds (2001), 
mainly due to the lack of low-metallicity galaxies observed 
at submm wavelengths because of the difficulty to perform 
such observations. Submm observations from Planck and newly 
available ground based submm instruments (e.g. SCUBA-2, 
LAB OCA) will be necessary to observe more low-metallicity 
environments at better sensitivity. Furthermore, the Herschel 
Guarantee Time Key Program, Dwarf Galaxy Survey (PI: S. 
Madden), is dedicated to the observation of a wide range of low- 
metallicity environments with the far-IR to submm PACS and 
SPIRE instruments. These observations will enable us to extend 
our knowledge on the D/G vs metallicity relation to lower metal- 
licities, using 70 to 500 /im data to constrain the peak and the 
submm slope of the SEDs. 



the quantity [Oii]^3727+[Oiii]i4959,15007 / US, also called the R; 
parameter (Pagel et al., 1979). 



4.2.3. Outlying galaxies 

In this §, we focus on the galaxies that deviate from the relation 
between the D/G and oxygen abundance. 

MCG+02-04-025 - The II2 mass could not be estimated for 
the galaxy due to the lack of CO observations for this galaxy. 
The estimated D/G is thus an upper Umit. 

M51 - This galaxy has an unusually "cold" global SED 
(Fig. 8). Applying only IRAS data, the SED modelling gives 
much lower total mass of dust but the use of Spitzer data, 
especially the 160 ixm of MIPS already significantly increases 
the dust mass, thus increases the D/G. There are still difficulties 
in estimating an average metallicity due to the extension of the 
galaxy and the fact that metallicity usually varies as a function 
of radius for spiral galaxies (e.g. PhilUpps & Edmunds, 1991; 
Issa et al., 1990). For M51, Moustakas & Kennicutt (2006) 
also estimated the metallicity to be 12+log(0/H)=8.7 using 
Pilyugin & Thuan (2005) but their estimate reaches 9.2 using 
the McGaugh (1991) strong-line calibration of the R23. A higher 
metallicity would, in fact, make the galaxy move rightward, 
close to the James et al. (2002) relationship but we prefer to keep 
the value 8.7 estimated by Bresolin et al. (2004), which is based 
on high-quality Hn region observations and calculated in the 
arms of the galaxy, where star-forming regions are concentrated. 
Nevertheless, M51 could have a strong gradient of metallicity 
that affects its position on this plot. 

M83 - The SED of this galaxy is plotted in Fig. 8. This 
galaxy is surrounded by a very large Hi envelope, with a 
radius > 125 kpc - 6.5 times the optical radius (Huchtmeier 
& Bohnenstengel, 1981). The total Hi mass of this envelope 
was estimated to be 2 x lO^" Mq (Tilanus & Allen, 1993). 
The Hi mass of Table 1 only accounts for that included in the 
photometric apertures used to estimate the fluxes we use in the 
SED modelling to derive the total mass of dust. The low D/G 
we derive is thus not due to an overestimation of the Hi mass. 
Lundgren et al. (2004) studied the kinematics of molecular gas 
of this galaxy with CO observations performed at the SEST 
and estimate the total gas mass, including H2 and Hi to be ~ 
7.8 X lO'^Mo, close to that used by Galliano et al. (2008) for 
the D/G estimate plotted in Fig. 7. The metalhcity used in this 
paper (12-i-log(0/H)=9.2) was determined by Webster & Smith 
(1983) using the R23 method used in our study. Engelbracht 
et al. (2008) recently re-estimated the oxygen abundance of 
M83 with the direct electron-temperature (Te) method applied 
to individual H n regions (Skillman et al., 1998) and found 
12-i-log(0/H)=8.62. This metallicity would thus bring the 
galaxy closer to the model predictions. 

IC342 - The submm constraint we used is the 1mm flux 
and is an upper limit (Fig. 8), meaning that the dust mass we 
estimated is also an upper limit too. The location of IC342 at 
low Galactic latitudes makes it difficult to study because of its 
large extent. It is affected by foreground extinction due to our 
Galaxy and the backgroimd subtraction issues it generates. This 
galaxy possesses a very peculiar optical/near-infrared shape 
(Galliano et al., 2008). 

He 2-10 - This Blue Compact Dwarf (BCD) differs from 
other BCDs because of its near solar metallicity (Kobulnicky 
et al., 1999) and contains a relatively large amount of molecular 
gas (Meier et al., 2001), which is unusual for low-metallicity 
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Fig. 7. a) Dust-to-gas mass ratio as a function of oxygen abundance for our restricted sample of galaxies that were observed at submm wavelengths. 
For the convention on symbols, refer to Fig. 4. Red symbols show the galaxies that present a submm excess and for which a cold dust (lOK) 
component with an emissivity fj=l is added. In this plot, the dust masses of the galaxies are derived from our SED models without introducing 
submm observational constraints. The solid line and the dotted line are the predictions of the dust formation models from Edmunds (2001). The 
gray stripe is the expectation of the one-zone, single-phase chemical evolution model of Dwek (1998) presented in Galliano et al. (2008). Both 
models are presented in § 2. The dashed line is the linear regression of the broad sample of §2 where all galaxies, including those without submm 
data, are used, b) Same but the dust masses of the galaxies are derived from our SED models using submm observational data. Labeled galaxies 
are discussed in § 4. 1 . Error bars refer to errors in the dust mass estimates and do not take the gas mass errors into account, c) We gather the two 
previous plots to emphasize how the D/G "moves" (dot-dash vertical lines) in the diagram when submm data are used to estimate the dust mass in 
the SED modelling. 
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Wavelength (^m) Wavelength (/im) 

Fig. 8. SED models of the galaxies M51, M83, IC342 and NGC 3190. The SEDs modelled with submm data are indicated by the black lines. 
Dashed lines indicates the model obtained when submm data are omitted. Observational constraints are superimposed (filled circles). When the 
error bars are not shown, the errors are smaller than symbols. The thick orange line shows the IRS spectrum used in the SED modelling when 
available. The open circles indicate the expected modelled fluxes integrated over the instrumental bands. The green and red lines respectively 
distinguish the stellar and dust contributions. 



dwarf galaxies. Nevertheless, the metallicity of the galaxy is 
rather uncertain due to the lack of measured [O III] /14363. 

NGC 3190 - No H2 mass estimates have been obtained in 
NGC 3190. The value presented on the plot is thus an upper 
limit. An amount of molecular gas equivalent to that of the Hi 
mass would lead to a higher D/G. The SED of this galaxy is 
plotted in Fig. 8 

Ham 11 - Using Spitzer data, Engelbracht et al. (2008) 
note that Haro 1 1 possesses an intriguing SED, dominated by 
the mid-IR emssion and showing weak emission at 70 fim. A 
D/G of ~ 10"^ would be expected, using the Galliano et al. 
(2008) chemical model, for a galaxy with such metallicity 
(12-i-log(0/H)=7.9). We find an unusually high dust/Hi mass 
ratio < 2 X 10"', using the Hi mass upper limit of 10** given 
by Bergvall et al. (2000). This high value could suggest that 
the scenario of cold dust may be physically inappropriate for 
this galaxy if the totalgas reservoir is estimated correctly. In the 
particular case of Haro 11, Galametz et al. (2009) show that, 
even taking all the possible upper limits of gas into account in 
the total gas mass estimate (H2 derived from CO(l-O) observa- 
tions in Bergvall & Ostlin (2002) or from [CII] measurement in 
Bergvall et al. (2000), ionized gas mass estimate in Ostlin et al. 
(1999)), the galaxy is an outlier of the Schmidt law (Kennicutt, 
1998), relation which relates the star formation rate of a galaxy 
to its total IR luminosity. They derived a "missing gas mass" of 
~ 10'" Mq, a value consistent with that independantly derived 



using the dust mass. A large amount of gas is thus not 'detected' 
by current observations in Haro 1 1 . 

The molecular gas mass, commonly inferred from the CO 
tracer, is very difficult to estimate in low-metallicity galaxies. 
The CO to H2 conversion factor is usually calibrated with cold (» 
12K) Galactic clouds (solar metallicity) and may not necessar- 
ily resemble the physical conditions of low-metallicity objects. 
Some studies show that this conversion factor can be far higher 
in low-metallicity environments than in metal-rich ones (Wilson, 
1995; Madden et al., 1997; Leroy et al., 2005, 2007). Due to 
the high excitation and the paucity of the dust in low-metallicity 
environments, CO is generally a poor tracer of molecular gas. 
The opacity of the H2 and CO lines can create a condition 
where the molecules are self-shielded, an effect which can be 
efficient for hiding H2 in regions where CO is photodissoci- 
ated. Cold dust could reside in small shielded molecular/cold 
dust cores embedded in [CII] emitting envelopes where UV ra- 
diation would hardly penetrate. High values of L([CII])/L(CO) 
have indeed been found in low-metallicity galaxies sometimes 
up to 10 times higher than normal dusty starbursts or spirals 
(PogUtsch et al., 1995; Madden et al., 1997; Madden, 2000; 
Hunter et al., 2001; Galliano et al., 2003). This high ratio of ob- 
served L([CII])/L(CO) generally suggests a clumpy ISM. It is, 
for instance, the case in the galaxy Haro 11 discussed above, 
with L([CII])A.(CO)> 10^ (Bergvall et al., 2000). CO is also 
known to be a poor tracer of dense gas. Further observations 
of the far-IR fine structure lines will be done with the Herschel 
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Dwarf Galaxy Survey (PI: S. Madden) to obtain a more accurate 
estimate of the gas content of dwarf galaxies. 



4.3. Conditions for submm excess 

We find that for 9 galaxies of our sample, a submm excess is 
detected (Fig. 3). To study the physical conditions leading to an 
excess emission at submm wavelengths, we plot the 850 fim flux 
as a function of the 160 fim flux (Fig. 9a). When the 160 fim or 
850 ijm fluxes are not directly observed, values are derived from 
the SED models we generated. Red symbols indicate galaxies 
where a submm excess is detected (SEDs modelled with an ad- 
ditional very cold dust component). While other galaxies follow 
a tight correlation, galaxies with submm excess show a 160-850 
correlation that seems to be shifted. Figure 9b shows the 160 
/vm / 850 fim flux density ratios as a function of metallicity. The 
160 fim to 850 fim flux ratio is lower for galaxies in which a 
submm excess is detected (most of the time below ~ 65, see the 
red symbols) but no clear correlation is found at this point with 
metallicity. Nevertheless, the sample is clearly lacking very low- 
metallicity galaxies (12-i-log(0/H)<7.9), which prevents us from 
studying the submm excess with an unbiaised metallicity cover- 
age. 

We also note that a submm excess is detected in dustier 
galaxies (NGC 4303 and NGC 2903). hi NGC 4303, the lack of 
far-IR measures between 100 and 450 jim does not enable us to 
correctly sample the peak of the dust SED and the submm slope. 
Observations sampling this wavelength range (with Herschel for 
instance) would be required to confirm the excess in this galax- 
ies. The dwarf galaxy He 2-10, which has a nearly solar metal- 
licity, also shows a submm excess. We note that our fit has dif- 
ficulties to fit the 100 and 160 jdm fluxes for this galaxy as weU 
as for NGC 2903. However, the 850 yum flux clearly exhibits an 
excess compared to the submm slope sampled by the 160 and 
450 jum fluxes. 

Submm excess is most evidently detected in low-metallicity 
galaxies. If we assume that this excess is due to very cold dust, 
could it be possible that very cold dust is more apparent in dwarf 
galaxies compared to more metal rich galaxies? Indeed, on av- 
erage, global SEDs of blue compact galaxies generally peak at 
relatively short wavelengths (30-60 jum), compared to more qui- 
escent, higher metallicity sources. For a metal-rich galaxy, peak- 
ing around 100 yum for example, the submm emission at 870 
jum should be very difficult to distinguish from the emission of 
warm grains while in the case of a low-metallicity galaxy peak- 
ing shortward, the separation between two types of grain popula- 
tion temperatures, one warm and one very cold, may be easier to 
detect. We test this idea by modelling four galaxies with differ- 
ent metalHcities (M51, M82, SBS 0335-052, Haro 11) with an 
independent very cold dust component, supposing that the mass 
and the properties (yS, T) of the very cold dust (VCD) is the same, 
thus Mvc/)/M^„s,=constant. These tests lead to the inverse con- 
clusion. The difference between the temperature of dust grains 
and the very cold grains is smaller in galaxies having higher 
metalHcities. With a constant mass ratio, the emissivity ratio be- 
tween the two components is thus higher in the dustier galaxies 
than in low-metallicity galaxies. If a very cold dust component 
exists with the same characteristics, our test here implies that 
it should be more evident as submm excess in dustier galaxies 
rather than in low-metallicity galaxies. 

The SEDs of low-metallicity galaxies usually show rela- 
tively "warm" (fvo/fieo > 1) dust temperatures. For these "warm" 
galaxies, the temperature of the SED is then well constrained but 
the flux at 870 fim will significantly influence the total amount 
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Fig. 9. a) Correlation between the 160 yum and the 850 //m fluxes of the 
galaxies. A linear regression is overlaid (dashed line). When the galax- 
ies do not possess an observational constraint at 850 yum, the flux value 
was deduced from the SED model, b) 160 yum / 850 yum flux density ra- 
tios as a function of metallicity. For convention on colors and symbols, 
see Fig. 4. Red symbols show the galaxies that present a submm excess 
and for which a cold dust (lOK) component with an emissivity y6=l is 
added. 



of dust while in "colder" galaxies (fvo/fieo < 1)> the SED model 
without submm data already leads to a large amount of dust. 
We plotted the ratio of the dust mass obtained with submm data 
over the dust mass without submm observations as a function of 
fvo/fieo but did not find any correlation for our galaxies. We also 
plot the wavelength of the IR peak as a function of metallicity in 
Fig. 10. We do not clearly detect a correlation between the posi- 
tion of the dust emission peak and the metallicity of the galaxy. 
However, we note that a submm excess is not detected when the 
SED peaks above 85 jjm. In some galaxies, the fact that a submm 
excess is not detected could also be due to an underestimate 
of the submm fluxes due to background subtraction issues (e.g. 
NGC 6946) or a lack of coverage of the complete emission in 
the submm observations (e.g. NGC 1097). Having several obser- 
vational constraints in the submm slope of the SED would also 
help us to be more confident in the presence of submm excess: as 
shown in Fig. 3, most of the galaxies where submm excess was 
detected were actually observed at different submm wavelengths 
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(at 450 and 850 yum). The addition of observational constraints 
will be a solution to avoid biases and sample a possible knee in 
the submm emission or a smoothness of the far-IR/submm slope, 
that could favor one of the two possible theories to explain the 
excess: a cold dust population or a change in the dust emissivity. 

Herschel is currently observing in the 70 to 500 yum wave- 
length range and will not only bring the lacking data to sample 
the SED but also the high resolution crucial to understand the 
origin of cold dust emission in galaxies. Herschel is also prob- 
ing the coldest phases of dust in high-z galaxies. Santini et al. 
(2010) recently found that the dust content of their high-z sub- 
millimeter galaxies appears to be far higher than that expected 
from their metallicity, probably revealing different dust proper- 
ties or dust growth mechanism in those environments. We cau- 
tion the accuracy of dust masses estimates made without rest- 
frame submm observations. Accurate dust mass estimates in 
high-z galaxies are necessary to enable us to study the evolution 
of metals through cosmic time. 
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Fig. 10. Wavelength of the peak of the dust emission as a function of 
metallicity. For convention on colors and symbols, see Fig. 4. Red sym- 
bols show the galaxies that present a submm excess and for which a 
cold dust (lOK) component with an emissivity yS=l is added. 



5. Discussion and Conclusions 

We have gathered observational constraints of a large sample 
of galaxies from low-metallicity to metal-rich environments to 
study the influence of submm measurements on the total dust 
mass estimates and the consequences on the D/G. We use mid- 
IR and far-IR data from the literature. Our second sample is con- 
fined to galaxies having submm data beyond 160 //m, namely 
850 or 870 yum, and sometimes 450 jum, leading to a broad sam- 
ple of 52 galaxies with submm observations. Submm observa- 
tions probe the coldest dust grain emission of galaxies, a com- 
ponent that could account for a significant fraction of the total 
dust mass of galaxies. 

We model the SEDs of our sample and derive their total dust 
masses with and without submm fluxes to determine the impor- 
tance of submm data in determining the dust SEDs. We find that 
submm data are crucial to accurately determine the dust masses. 
The error bars on the total dust mass estimates significantly de- 
crease when using submm constraints. We show that submm data 



systematically affect the dust mass derived for our galaxies for 
different reasons: a) for dustier galaxies for which the SED usu- 
ally peaks at longer wavelengths, submm fluxes are crucial to 
position the peak and the Rayleigh- Jeans slope of their SED and 
avoid an overestimation of the total dust mass of the galaxy usu- 
ally observed when the SED modelling is performed using data 
at wavelengths only as long as 160 /vm; b) the submm wave- 
length domain sometimes harbours an excess that may imply a 
large amount of very cold dust. The lack of submm constraints 
will lead, in those cases, to an underestimate of the total dust 
mass of the galaxies. We caution the estimate of dust masses 
made without submm observations. 

We find that for 9 galaxies of our sample, mostly low- 
metallicity galaxies, our fiducial model can not fit the mid-IR 
to submm observational constraints and a submm excess is de- 
tected. We model this excess with a very cold dust component, 
using a temperature of lOK and a dust grain emissivity index 
of 1 . The early results of Herschel also highlight a submm ex- 
cess in dwarf galaxies (O'Halloran et al., 2010; Grossi et al., 
2010). First results based on observations with the Planck mis- 
sion in nearby galaxies have also lead to evidence of cold dust 
(Planck Collaboration, 201 Ic). In the Magellanic Clouds, Planck 
Collaboration (201 lb) have lead to the detection of an excess at 
submm/mm wavelengths. Average emissivity indexes /? of 1.5 
and 1 .2 were found for the LMC and SMC respectively, so flat- 
ter than observed in the Milky Way. Even taking the flatter slope 
into account, there is still evidence for a further flattening in the 
submm regime for the SMC, that has been modeled with a com- 
bination of thermal and spinning dust emission. 

Gathering the gas mass information for our galaxies, we es- 
timate the D/G of our sample. We find a tightened relation be- 
tween the D/G and the metallicity of the galaxies when the mod- 
els take into account submm observations. This effect is present 
in spite of our model assumption of cold dust for the submm 
excess. 

More submm observations of low-metallicity galaxies, 
which will be obtained with Herschel and Planck will extend the 
coverage to a broader range of metallicity and enable us to study 
how the D/G evolves in more metal poor environements. The 
spatial resolution and sensitivity of Herschel combined with cur- 
rent ground based instruments (for example, JCMT / SCUBA-2; 
APEX/Laboca) and ALMA will help us probe the submm emis- 
sion of nearby galaxies in detail to better isolate the excess emis- 
sion and understand the physics of the coldest phases of dust 
and the chemical enrichment thus drawing the global picture of 
the dust distribution within galaxies. While presently the nature 
of this submm excess seen in low metallicity environments is 
still not understood, the higher spatial resolution submm obser- 
vations will help us disentangle the different scenarios: emission 
from very cold dust, possible changes in the dust grain emissiv- 
ity, spinning dust etc. 
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